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STABILITY  OF  LAMINAR  ELECTRON  LAYERS 

I.  Introduction 

There  exists  by  now  a  truly  enormous  literature  on  the  subjects  of  the 
equilibrium  and  stability  of  layers  of  charged  particles  In  various 
geometries.  The  earliest  studies*-^  were  conducted  in  order  to  understand  the 
"slipping  stream"  or  dlocotron  effect  In  the  operation  of  the  first 
magnetrons.  Somewhat  later  the  importance  of  curvature  effects  was  realized 
when  beams  became  relativistic,  as  In  particle  accelerators;  the  resulting 

o  q 

"negative  mass"  instability  ’  completely  dominates  the  planar  beam  dlocotron 

effect  at  sufficiently  high  energies  (only  a  few  tens  of  keV  in  many  practical 

devices  of  interest).  The  negative  mass  effect  was  also  investigated  in 

connection  with  some  controlled  fusion  research  devices and  other 
1 3 

machines . 

In  recent  times  there  has  been  renewed  interest  in  high  power,  high 
efficiency  microwave  devices  as  well  as  in  accelerators  capable  of  high 
current  operation.  Spurred  by  the  discovery  of  the  electron  cyclotron  maser 
(or  gyrotron)  effect, research  in  the  field  of  short  wavelength,  high 
power  microwave  devices  requiring  no  slow  wave  structure  has  been  vigorously 
pursued.*^  Operation  of  these  new  devices  depends  fundamentally  on  the 
negative  mass  effect  as  enhanced  by  a  synchronism  of  the  particles'  angular 
motion  with  the  temporal  and  angular  variation  of  a  "cold"  waveguide  mode. 

This  enhancement,  though  not  called  the  maser  effect,  was  first  noted  in  the 
classic  work  of  Briggs  and  Neil*®  and  has  been  further  elaborated  in  ref.  18. 

The  acceleration  of  large  currents  of  electrons  is  a  formidable  problem 
which  has  also  received  considerable  attention  recently.  In  cyclic  devices  it 
is  possible  to  construct  high  current  equilibria  in  many  cases  but 

various  instabilities,  the  negative  mass  instability  prominent  among  them,  may 
limit  achievable  currents  to  smaller  values  than  equilibrium  considerations 
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alone  would  suggest.  In  these  devices  the  effects  of  self  fields  on  the 
negative  mass  Instability  and  on  the  stabilization  mechanism  become  Important 
to  consider;  In  one  device  these  effects  have  been  predicted  to  lead  to  a 
peculiar  double-valued  feature  In  the  current  vs.  energy  spread  stability 
curve.  It  was  this  particular  result  which  Initially  prompted  the  study  of 
self  field  effects  reported  here  for  a  simpler  (2D)  model. 

The  model  we  consider  consists  of  a  layer  of  charged  particles  (we  will 

think  of  them  as  electrons  but  Ions  may  be  trivially  substituted)  moving  in 

circles  about  a  common  axis,  as  shown  in  Fig.  1.  This  restriction  to  laminar 

or  "cold"  electron  flow  will  tend  to  overestimate  actual  growth  rates  of  the 

modes  we  study  since  the  effects  of  betatron  oscillations  and  of  nonaxis- 

encircling  particles  are  stabilizing.  An  analysis  including  these  effects  is 

properly  done  in  phase  space  using  the  correct  equilibrium  distribution 

function;  such  an  analysis  is  significantly  more  complex  than  that  given  here 

and  may  be  carried  out  only  in  an  approximate  way.  The  laminar  flow  case 

includes  all  essential  physics  and  has  great  simplicity,  allowing  an  exact 

treatment  of  the  linearized  problem,  to  recommend  it.  The  equilibrium, 

discussed  in  the  following  section,  is  supported  by  a  combination  of  self  and 

externally  applied  radial  electric  and  axial  magnetic  fields.  We  impose  no 

a  priori  restriction  on  the  relative  magnitudes  of  the  three  terms — 

centrifugal  force,  electric  force,  v  x  S  force — in  the  equilibrium  force 

23 

balance.  This  configuration  is  thus  a  reasonable  model  for  the  Astron, 
magnetron,  gyrotron,  orbitron,^  peniotron^  and  heliotron^  and  includes 
correctly  the  crucial  radial  and  azimuthal  particle  dynamics  found  in  particle 
accelerators.  This  "E-layer"  model  has  the  virtue  that  the  linear  stability 
problem  may  be  formulated  exactly,  for  arbitrary  particle  energies,  including 
all  effects  of  self  fields  and  all  effects  of  relativity  and 


2 


electromagnetism.  It  is  perhaps  a  bit  surprising  that,  despite  the 

venerability  of  the  topic,  this  exact  formulation  has  not  been  carried  out 
27 

earlier.  The  eigenvalue  problem  governing  the  stability  of  the  E-layer  is 

derived  and  analyzed  in  Section  III,  below. 

The  desirability  of  a  completely  general  treatment  encompassing  many 

devices  and  a  large  parameter  space  is  related  to  the  ease  with  which  various 

familiar  results  for  special  cases  may  be  recovered  in  appropriate  limits.  A 

dispersion  relation  for  the  so-called  longitudinal  mode  obtained  analytically 

in  the  thin  beam  limit  from  the  ordinary  differential  equation  governing  the 

RF  field  reproduces  all  standard  special  case  dispersion  relations  (negative 

mass,  electron  cyclotron  maser,  diocotron)  in  a  straightforward  way;  in  doing 

so,  a  puzzle  is  resolved  concerning  the  survival  of  a  finite  negative  mass 

growth  rate  in  the  planar  limit  and  a  method  suggests  itself  on  how  either  to 

28 

maximize  or  eliminate  altogether  the  negative  mass  instability  growth. 
Extensive  testing  of  the  dispersion  relationship  against  a  numerical  solution 
of  the  eigenvalue  problem  shows  excellent  agreement,  as  reported  in  Section 
III.  The  classical  diocotron  dispersion  relation  is  discussed  and  recovered 
from  our  formalism  in  Section  IV. 

Our  treatment  of  the  transverse  mode,  which  has  been  invoked  in  the 
theory  of  the  penlotron  device  is  somewhat  less  comprehensive.  This  mode  is 
unstable  only  when  electron  motion  is  synchronous  with  a  cavity  mode  and  so  is 
of  a  more  specialized  nature  than  the  negative  mass  or  diocotron  instabilities 
which  arise  from  intrinsic  properties  of  geometry  and  shear  flow,  not  from 
interaction  with  external  structures.  Still,  the  utility  of  the  transverse 
mode  for  microwave  generation  may  not  yet  have  been  fully  exploited.  Ve 
comment  on  some  features  of  this  mode  in  Section  V. 


II.  Cold  E-Layer  Equilibrium 

We  consider  an  idealized  model  of  an  electron  layer  in  which  all 
electrons  circle  a  common  axis  as  shown  in  Fig.  1.  We  will  neglect  the 
effects  of  betatron  oscillations  and  of  axial  motion  in  both  perturbed  and 
unperturbed  states.  Inclusion  of  betatron  oscillations  or  axial  perturbations 
is  expected  to  have  a  stabilizing  influence  on  the  collective  modes  we  will 
study.  We  further  assume  that  in  general  the  layer  is  enclosed  in  a  coaxial 
waveguide  with  smooth,  perfectly  conducting  walls  at  r*a  and  r*b,  as  shown  in 
the  figure.  Absence  of  either  or  both  walls  results  only  in  a  change  of 
certain  boundary  conditions. 

Though  this  model  is  simple  to  describe  it  is  surprisingly  rich  in 
content.  Depending  on  beam  and  geometric  parameters  and  types  of  applied 

OO  |  7 

fields,  it  may  be  taken  as  a  good  description  of  the  Astron,  gyrotron, 

24  2S  2fi 

orbitron,  peniotron,  heliotron  °  and  cross  field  microwave  devices  and  may 

o  a 

also  be  of  some  interest  in  accelerator  '  and  space  physics  stability 
problems  as  well  as  in  the  theory  of  magnetic  insulation. 

A  A 

The  equations  governing  the  equilibrium  quantities  vo(r)9,  EQ(r)r, 

A 

BQ(r)z  for  a  specified  density  profile  nQ(r)  are 


y  v 
'o  o 


-  [E  +  v  B  ] 
m  1  o  o  oJ 


(1) 
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r  dr 
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-en  /e 
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where  -e  and  m  are  the  electron  charge  and  mass,  y  is  the  usual  relativistic 
2  -  1/2 

factor,  (1-(vq/c)  )  ,  c  is  the  speed  of  light  in  vacuum,  and  eQ  and  are 
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Che  permittivity  and  permeability  of  free  space.  In  addition  to  the  density 
profile,  we  take  as  given  the  total  electrostatic  potential  difference  and  the 
total  magnetic  flux  contained  between  r*a  and  r-b. 

Using  these  conditions  the  solution  to  Eq.  (2)  is  immediate  and  Eqs.  (1) 
and  (3)  may  be  combined  to  yield  a  single  differential  equation  for 

U  5  6oYo  S  voVc: 


; 

u 


1  ~  5/y2  ~  y2h 
u  r  'o  lo 

r  1  1  +  h 


(4) 


2  z  2  2 

where  £  ■  w  /&>  ,  us  ■  e  n  /my  e  is  the  plasma  frequency  of  the  layer, 
pop  O  0  0 

uso  ■  v0/r»  and  a  prime  denotes  d/dr.  In  Eq.  (4)  we  introduce  the  important 
quantity  h,  defined  by, 


h 


erE 


2  2  3 

mc  80  ’'o 


(5) 
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which  is  (l/yQ)  times  the  ratio  of  the  electric  force  to  the  centrifugal  force 
experienced  by  an  electron  in  equilibrium  at  radius  r. 

Equation  (4)  must  be  solved  subject  to  an  initial  condition  u(r^) 
satisfying  Eq.  (1)  which  we  rewrite  in  the  dimensionless  form 


u2  -  o^+l)172  +  aBu 


(6) 


where  a-  ■  erE  /me  end 
E  o 


erB  /me. 
o 


Og  is  known  at  r^  but  ag  is  not.  The 


entire  equilibrium  problem  then  reduces  to  choosing  otg(r^)  such  that  the 


resulting  total  flux 


b 

*  -  2irf  dr  rB  (7) 

O'  o 
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i9  a  specified  number.  A  numerical  method  for  doing  this  is  described  in 
Appendix  A.  Here  we  briefly  consider  the  character  of  some  of  the 


solutions.  The  equilibrium  thus  constructed  will  be  used  in  the  numerical 
solution  of  the  eigenvalue  problem  in  Section  III. 

Given  c^,(r 

solutions  for  u(r,),  depending  on  the  values  of  a-,  and  a„;  without  loss  of 
generality  we  take  ^(r^)  >  0-  The  situation  in  the  a£-  ag  plane  is  then 
depicted  in  Fig.  2.  By  a  positive  root,  denoted  by  (+)  in  the  figure,  we  mean 
a  right-handed  rotation  about  $  (clockwise,  looking  in  the  direction  of  1$) . 

The  equation  for  the  boundary  curve  C  in  Fig.  2  is  the  condition  for  Eq.  (6) 
to  have  a  double  root: 


and  ag(r^),  Eq.  (6)  may  have  0,  1,  or  2  real  valued 


“E 


*  * 
u  (u 


(U*2+  1)1/2 


(8) 


where  u*  -  (8/3)1/2sinh[-j  sinh~1[^(3/2)3/2aB] ] .  On  C,  1  +  h  -  0. 

There  are  three  basic  types  of  equilibria  found  in  devices  of  practical 
interest:  Type  I,  in  which  a  magnetic  field  is  used  to  balance  centrifugal 
force  (<*E  is  small,  e.g.  Astron,  gyrotron,  particle  accelerators);  Type  II,  in 
which  electric  and  magnetic  forces  balance  (inertia  is  small,  e.g.  crossed 
field  microwave  devices);  and  Type  III,  in  which  an  electric  field  is  used  to 
balance  centrifugal  force  (ag  is  small,  e.g.  orbitron,  heliotron).  Type  I  may 
have  only  "+"  roots  while  Types  II  and  III  may  have  either  "+"  or  roots. 

In  the  following  analysis,  however,  we  make  no^  assumption  regarding  the 
relative  magnitudes  of  the  electric  field,  magnetic  field,  and  inertia  terms. 

It  is  reasonable  to  ask  whether  stability  considerations  favor  one  type 
of  equilibrium  over  another  for  a  layer  of  given  kinetic  energy.  To  our 
knowledge  there  has  been  no  definitive  answer  to  thi9  fundamental  question. 
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For  radiation  source  applications  one  arguably  might  want  the  most  unstable 
configuration  while  for  other  applications,  like  accelerators,  one  would  want 
the  most  stable  one.  We  proceed  to  examine  the  stability  of  a  general 
equilibrium  describable  by  Eqs.  (1)  -  (3)  in  order  to  address  this  question, 
our  goal  being  a  complete  parametric  study  of  the  small  signal  behavior. 


III.  Stability  of  the  Layer 

We  consider  perturbations  on  the  equilibrium  described  in  the  previous 
section.  The  equations  governing  the  layer  are,  simply 


^ftT  +  m 

-  -(£  +  vx$) 
m 

(9) 

$x£  = 

_  3$ 
at 

=  -en/e 

0 

(10a, b) 

- 

l  a£  •> 

~2  at  '  Vnv 

=  0 

(11a, b) 

where  n  is  the  number  density. 

Our  discussion  will  be  limited  to  consideration  of  perturbations  in 

A  A 

the  r  and  0  directions,  for  reasons  described  below.  Writing  all 
quantities  ?  as 

nr.O.t)  -  ?o(r)  +  ei(',,t_i0)'?1(S')(r;w)  (12) 

i l 

where  fQ(r)  is  the  equilibrium  value,  and  retaining  only  terms  linear 

(l) 

in  one  finds  that  the  linearized  versions  of  Eqs.  (9)  -  (11)  decouple 

into  two  sets  governing  {vlr »vl0,Elr ,El0,Blz .n^  and  { vlz -Elz » Blr » Bl0}  which 
we  identify  as  TE  modes  and  TM  modes  respectively.  (Here  and  below  we 
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•  4 


•  4 


•  4 


•  4 


•  < 


(i) 

write  (r;w)  simply  as  for  all  first  order  quantities.)  The  TM  modes 

may  easily  be  shown  to  be  neutrally  stable  for  perfect  conductor  boundary 
conditions  (damped,  for  resistive  wall  boundary  conditions).  They  represent 
simple  oscillations  of  the  electrons  along  the  z  direction  in  response  to  the 
cavity  mode  fields;  the  equilibrium  model  does  not  provide  any  free  energy  to 
excite  the  TM  modes.  The  TE  modes,  on  the  other  hand,  are  potentially 
unstable.  If  we  had  allowed  a  finite  axial  wavelength  for  the  perturbed 
quantities  the  TM  and  TE  modes  would  be  coupled,  but  it  has  been  found10  that 
the  effects  of  finite  axial  wavelength  are  stabilizing  for  the  coupled  case. 
Consequently  we  focus  attention  exclusively  on  the  TE  waves. 

The  Euler  equation  (9),  upon  linearization,  gives  the  fluid  response  in 
terms  of  the  RF  fields  of  the  TE  wave: 


Vie 


—  [E,  +  v  B,  1 
my  1  Ir  o  lzJ 
o 


(13a) 


“o^ir^le 


(13b) 


where  we  have  defined  H  =■  ur-£oi  ,  P  =  y  (1+h),  Q  *  v  /u>  +  h,  and 

o  o  00 

v'=  ^  vq.  Using  the  solution  of  Eqs.  (13a, b)  in  the  definitions  of  the  RF 
currents  and  and  combining  the  linearized  Maxwell's  equations  (10)  - 
(11)  allows  one  to  write  a  single  differential  equation  for  4>  =  rE^: 


4<tXJF>  + ■ 0 


where 
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(15) 
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To  0 


2  2 
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(17) 


_  .2,  2  ,  2  ,  2,2  2 ,  2 
D  *  PQ  -  JJ  /to  =  l-E/y  +  Y  h  -  H  /id 
o  oo  o 


(18) 


and  other  symbols  have  been  defined  previously.  [Equation  (14)  is  identical 
to  Eq.  (4)  of  Ref.  28,  except  for  a  few  differences  in  notation.] 

The  other  RF  fields  are  given  in  terms  of  $  by 

1  .  i  *\  tlKJ-  1 

Elr  *  lrA  L<7  +  5D“l0v  /O*' - 1  (-D  (1-rti)*]  (19) 

c 

Blz  -  -  ^  A~1[(o)  +  ^D"1^)^  -  |  o>oCD'1(1-H1)^]  (20) 

c 

and  the  perturbed  velocities  may  be  easily  recovered  using  Eqs.  (19)  and  (20) 
in  Eqs.  (13a, b);  the  perturbed  density  is  obtained  from  n^>*-(co/e)?»S^. 

Equation  (14)  must  be  solved  subject  to  $(a)  *  d> ( b)  ■  0  (if  the  walls  are 
perfect  conductors)  in  order  to  obtain  the  eigenvalue,  o>.  Hereafter  we  shall 
assume  that  the  electron  layer  is  restricted  to  the  annular  region  r^<r<r2 
[Fig.  1].  In  this  case  it  is  convenient  to  formulate  the  problem  using  the 
explicit  vacuum  solutions  at  r»r^  and  r*^  anc*  t0  matc^  the  appropriate 
logarithmic  derivatives,  as  is  often  done  in  microwave  tube  theory,  in  order 
to  isolate  the  effect  of  wall  boundary  conditions.  Defining,  then,  the 


normalized  admittances  evaluated  just  outside  of  the  layer,  we  have 
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We  stress  that  the  normalized  wave  admittances  b+  and  b_  depend  only  on  the 
waveguide  geometry  exterior  to  the  electron  beam.  They  are  independent  of  the 
beam  or  its  dynamics.  b+  and  b_  are  evaluated  explicitly  in  Appendix  B  for 
some  practical  cases  of  interest. 

Let  us  examine  Eq.  (14)  more  closely.  In  the  complex  r  plane  Eq.  (14) 
has  singularities  at  points  where 


r-r. 


(21b) 


o)  -  £<jl)  -  0 

o 

(22) 

D  +  5/y^  -  0 

(23) 

A  ■  0. 

(24) 

The  first  of  these  clearly  represents  a  match  between  the  mode  frequency  and  a 
harmonic  of  the  particle  "cyclotron"  frequency.  Such  a  match  is  present  in 
the  negative  mass  instability,  cyclotron  maser  instability,  and  diocotron 
instability  all  of  which  are  described  by  Eq.  (14).  Indeed,  these  may  all  be 
considered  to  be  the  same  instability  in  this  sense, though  the  individual 
names  are  still  useful.  The  negative  mass  instability  is  fundamentally  a 
rotational  effect.  The  classical  explanation  attributes  it  to  the  decrease  of 
circulation  frequency  with  particle  energy  leading  to  growth  of  azimuthal 
bunches.  The  negative  mass  instability  operates  without  regard  for  the  cavity 


10 


modes  of  the  vacuum  chamber;  the  fields  are  not  well  approximated  by  the 
cavity  mode.  The  instability  is  strongly  enhanced  however  if  Re(<o)  happens  to 
be  a  cavity  mode,  as  first  pointed  out  by  Briggs  and  Neil.*®  In  recent  years 
this  synchronous  case  has  been  given  its  own  name,  the  cyclotron  maser 
instability,  mainly  in  the  literature  of  gyrotron  research.  » 10 

The  diocotron  instability,  originally  studied  in  connection  with  the 
development  of  the  first  cross  field  microwave  devices,  is  a  "residual" 
instability  in  this  context.  It  is  what  survives  in  the  non-relativistic 
and/or  planar  limits.  The  perturbed  electric  fields  are  basically 
electrostatic  in  nature  and  are  strongly  localized  at  the  position  of  the 
layer.  The  equation  governing  the  diocotron  instability* »^»6.7»29  may 
recovered  formally  from  Eq.  (14)  by  taking  c-m».  It  is  discussed  further  in 
Section  IV,  below. 

A  mode  satisfying  Eq.  (22)  is  sometimes  called  the  longitudinal  mode 
since  the  major  effect  on  the  beam  is  azimuthal  bunching.  Equation  (23),  on 
the  other  hand,  describes  what  has  been  called  the  transverse  mode  which 
involves  little  bunching,  but  significant  transverse  (radial)  motion  of  the 
beam.  If  h  is  small,  Eq.  (23)  becomes 

u  -  £u)q  ■  ±  a)Q*  (25) 

This  transverse  mode  has  been  invoked  to  explain  the  operation  of  the 
25 

peniotron  J  when  ai  corresponds  to  a  cavity  mode.  In  the  nonsynchronous  case 
the  transverse  mode  is  stable  in  the  absence  of  resistive  walls,  as  pointed 
out  in  ref.  10. 

Finally  we  note  that  Eq.  (24)  may  be  loosely  associated  with  an 
electromagnetic  mode.  In  vacuum  only  (24)  survives  as  a  singularity  but  in 


the  vacuum 


fact  it  may  be  shown  to  be  only  an  "apparent”  singularity;^0 
eigenfunctions  are  analytic  in  r,  except  possibly  at  r*0. 

We  proceed  to  analyze  Eq.  (14)  to  uncover  the  parametric  dependences  of 
the  growth  rates  of  the  unstable  modes.  In  order  to  make  progress 
analytically  we  consider  the  case  of  a  beam  of  uniform  density  and  of 
thickness  x  =  ~  r^  which  is  much  less  than  its  average  radius 

R  =  y  +  ri)*  ^  the  beam  is  sufficiently  thin  it  is  possible  to  consider 
a  Taylor  series  solution  about  r«R  of  Eq.  (14),  taking  care  to  check  at  the 
end  that  the  singularities  are  sufficiently  far  away  in  the  complex  r  plane, 

•J1 

so  as  not  to  disturb  the  series  convergence;  carrying  out  this  program  we 
find  the  dispersion  relation 

C  1  A1  1 
G+  -  G_  +  +  £  (G++  G_)(~  +  ±)  +  G+Gj 

R  Aq  o 

-  ,  C 

-  |({r<G+  -  GJ(^)  -  0  (26) 

o 

2 

correct  to  order  (x/R)  .  In  Eq.  (26)  G_,  G+  are  the  radial  logarithmic 
derivatives  of  <)>  evaluated  just  inside  the  layer  at  r*r^,  r**^  respectively, 
obtained  by  integrating  Eq.  (14)  across  the  beam  edges: 


- ^ -  [*  K  +  1^2  !>] 

t2,lA(r2,l) 


(27) 


where 


q(r) 


u  (n  i.  g(i  +  h) 

n  l,u>  ~  ad 


(28) 


and 
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(29) 


A(r)  3  Aq  +  (r-R)A1  + 

2(r)  *  C  +«•• 
o 


near  r»R. 


Expanding  G±  to  order  t/R,  Eq.  (26)  then  reduces  to 


•  4 


•  I 


|(b+  +  b _)  +  I  [-Rq'  -  q)(|b_+  q)  +  C  1-0  (30) 


•  4 


where  now  q  and  q'  are  understood  to  be  evaluated  at  r*R  and  where  in  Eq.  (30) 
we  have  (temporarily)  kept  terms  only  to  0(t/R).  No  assumption  has  been  made 
about  the  E-layer  density  or  current;  we  have  assumed  only  that  the  beam  is 


•  4 


"thin". 

Let  us  first  consider  the  longitudinal  mode,  that  is,  we  look  for  a  root 
2  2 

of  Eq.  (30)  with  In  |«u>o.  the  terms  within  the  brackets  in  Eq.  (30)  only 
/  2  -2 

the  q  and  q  terms  behave  as  for  small  Q.  Keeping  only  such  terms,  some 
algebra  gives  ^^=^(0): 

.  <■  iuj  (62  +  2h)  , 

al  5  TTTtV)  12)  +  <3l> 

+  -  1  +  y  h 

'o 

which,  for  h  ■  0,  is  the  classical  negative  mass  dispersion  relation  for  the 
Astron  configuration. 10  (The  0(*t/R)2  term  is  displayed  fully  in  Eq.  (36), 
below.)  We  remark  that  Eq.  (31),  including  the  second  order  terms,  can  be 
obtained  from  the  much  simplified  form  of  Eq.  (26), 


•  i 


•  4 


G+-  G_  +  tG+G_  -  0 


(26') 


where  terms  behaving  as  n  have  been  dropped  from  Eq.  (26). 


•  4 
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Ignoring  the  0(ft/R;  term  in  Eq.  (31)  for  a  moment,  we  may  interpret  the 

dispersion  relationship  (31)  as  follows.  As  we  have  seen,  the  factor 

(b+  +  b_)  represents  the  effects  of  the  container  structure,  and  the 
2 

factor  (t/R)(d^  is  proportional  to  the  beam  current.  These  two  terms  always 
appear  in  any  dispersion  relation  involving  an  electron  beam.  All  dynamical 
effects,  including  those  due  to  self  fields,  are  contained  in  the  factor 


(S*  +  2h) 

M  = - 2 - __  (32) 

(1  +  yV) 


for  a  sufficiently  thin  beam.  Some  insight  into  the  meaning  of  this  factor 

may  be  had  if  we  note  that  for  small  equilibrium  densities,  £,  it  is 

2 

proportional  to  do^/dc  where  e  ■  me  yq  “  e$0  18  the  total  energy  of  a  particle 
in  the  equilibrium: 


~  “08q  +  2h  *  Hit 

me2  1  +  y2h2-  £/y2 


(33) 


At  least  when  £  is  small  then  we  find  agreement  with  the  classical  explanation 

of  the  negative  mass  effect:  A  rotating  beam  is  unstable  (stable)  if  its 

equilibrium  rotation  frequency  is  a  decreasing  (increasing)  function  of  its 

(total)  energy,  that  is,  M  is  proportional  to  the  effective  azimuthal  Inertia 

(mass).  We  stress  that  for  finite  £  this  interpretation  begins  to  break  down; 

Eq.  (31)  was  derived  assuming  £  was  finite  yet  £  does  not  appear  in  the  factor 

M  but  clearly  does  appear  in  Eq.  (33). 

M,  though  simple  in  form,  has  many  interesting  features.  Let  us  consider 

its  dependence  on  h,  illustrated  in  Fig.  3.  Perhaps  the  most  interesting 

2 

property  of  M  is  that  it  experiences  a  change  in  sign  at  h  ■  -  BQ/2.  The 
possibility  for  using  this  property  of  equilibria  supported  by  a  radially 


14 


inward  electric  field  in  addition  to  a  magnetic  field  in  order  to  suppress  the 

OO 

negative  mass  instability  has  been  discussed  elsewhere.  We  note  that  since 

the  factor  M  is  independent  of  the  beam  current  or  of  the  waveguide  dimensions 

and  represents  a  purely  dynamical  quantity,  this  stabilization  of  the  negative 

mass  instability  by  a  radial  electric  field  is  expected  to  be  valid  even  in 

toroidal  geometry,  even  for  high  current,  very  cold  beams.  (The  method  is 

limited  to  use  for  moderately  low  energy  beams,  in  the  MeV  range  in  practice, 

since  the  applied  electric  field  required  to  change  the  sign  of  M  from 

3 

negative  to  positive  is  proportional  to  yo«  as  follows  from  the  definition  of 

h,  Eq.  (5).)  Note  that  there  is  a  "most  stable”  configuration  at  h«-l ,  at 

which  point  certain  singular  parts  of  C  in  Eq.  (14)  vanish  identically. 

Perhaps  equally  interesting  for  microwave  generation  applications,  is  the 

2 

occurrence  of  a  "most  unstable”  point  in  Fig.  3  at  h  ■  1/y  •  If  we  recall 

o 

2 

that  y  h  is  the  ratio  of  the  radial  electric  force  to  the  centrifugal  force  we 

2 

observe  that  the  choice  YQh  *  1  describes  a  configuration  in  which  the 

equilibrium  is  supported  solely  by  an  electric  field  (a  Type  III  equilibrium, 

in  the  language  of  Section  II).  This  result  suggests  that  for  a  given  beam 

energy,  a  microwave  source  such  as  the  orbitron^  in  which  an  annular  beam 

circles  a  positively  charged  wire  with  no  applied  magnetic  field,  might  have 

some  advantage  over  more  conventional  devices  like  magnetrons  (h<0),  inverted 
2 

magnetrons  (h»l/YQ),  or  gyrotrons(h»0) .  This  finding  is  novel.  We  remark 

2 

that  this  negative  peak  in  M  is  sharpest  at  low  energies  (small  SQ)- 

2  2 

In  the  non-relativistic  limit,  Bo+  0,  yo+  1,  the  sign  of  M  is  determined 
directly  by  the  sign  of  h,  which  remains  finite  in  this  limit.  Usually  we 
are  accustomed  to  thinking  that  the  negative  mass  instability  should  vanish  In 
the  non-relativistic  limit  when  there  is  no  gradient  in  the  magnetic  field. 
Here  we  see  however  that  an  equilibrium  electric  field  can  affect  the  sign  of 


Che  azimuthal  inertia  just  as  a  gradient  in  B  can.  If  we  define  an  equivalent 


field  index,  ne<j,  by 


T- 


(34) 


then 


n 


eq 


(Y^-Dh 


+  h)‘ 


(35) 


which  vanishes  when  h  does. 

One  final  comment  about  h  must  be  made.  The  electric  field  appearing  in 
the  definition  of  h,  Eq.  (5),  is  due  both  to  the  electric  charge  of  the  layer 
and  to  any  externally  applied  bias  potential.  When  the  contribution  of  the 
electric  charge  dominates  there  is  some  question  as  to  what  value  of  h  to  use 
in  Eq.  (31)  when  calculating  growth  rates  since  Eq  and  therefore  h  will  change 
sign  somewhere  within  the  layer.  For  this  reason,  in  all  numerical  examples 
considered  below  we  will  assume  that  h,  when  it  is  non-negligible ,  is 
dominated  by  the  contribution  of  an  externally  maintained  bias  potential. 

The  analytical  results  presented  up  to  this  point  are  subject  to  test  by 

numerical  integration  of  Eq.  (14)  subject  to  suitable  boundary  conditions.  We 

have  written  a  program  to  carry  out  this  task.  Given  a  density  profile,  the 

electrostatic  potential  difference  and  the  total  magnetic  flux  contained 

between  r«a  and  b,  the  program  calculates  all  equilibrium  quantities  then, 

given  the  mode  number  l,  locates  an  eigenvalue  u>.  In  all  examples  below  the 

density  profile  is  parabolic;  we  expect  good  agreement  with  Eq.  (31),  derived 

2  2  2 

using  a  flat  profile,  if  we  identify  (t/R)  uip  with  2(v/Yq)(c  /R  )  where  v  is 
Budker's  parameter,  the  number  of  particles  per  unit  length  times  the 


classical  electron  radius.  All  examples  use  a  perfectly  conducting  wall; 

effects  of  a  small  wall  resistivity,  which  may  be  very  important  under  some 

conditions,  are  discussed  briefly  toward  the  end  of  Appendix  B. 

Though  we  shall  eventually  test  Eq.  (31)  for  all  parameter  dependences, 

let  us  first  check  the  interesting  dependence  on  h,  discussed  above,  against  a 

numerical  solution  of  Eq.  (14).  Figure  4  presents  some  typical  results.  Here 

we  have  illustrated,  for  two  different  values  of  v/yq  a  comparison  of 

the  l  ■  1  mode  growth  rates  for  a  thin  (t/R  *.02)  beam  as  predicted  by  Eq. 

(31)  (dashed  line)  and  by  a  solution  of  Eq.  (14)  (solid  line).  (To  be 

precise,  what  is  plotted  on  the  numerical  solution  curve  is  Im(u))/uo(r^) . ) 

2 

The  transition  from  stable  to  unstable  behavior  occurs  at  h  ■  -  3Q/2 

(a  -  0.278  for  this  case)  Independent  of  £.  This  is  an  important  feature  of 

Eqs.  (31)  and  (32)  which  has  been  used  to  argue^®  that  the  stabilization 
2 

condition,  h<-BQ/2,  for  the  negative  mass  instability  is  independent  of  the 

beam  current.  Note  that  this  statement  cannot  be  made  on  the  basis  of  single 

particle  orbit  theory  alone  [cf.  Eq.  (33)];  the  stabilization  condition  is 

obtained  from  collective  mode  considerations  including  self  field  effects. 

Actually  a  small  amount  of  residual  growth  remains  for  h  less  than  but  very 
2 

close  to  -  3  /2,  a  feature  discussed  in  more  detail  in  the  following 
o 

section.  Agreement  between  Eq.  (31)  and  the  numerical  solution  is  best  at 

small  growth  rates  which  is  reasonable  if  we  recall  that  terms  behaving 

—1  ~2 
as  n  were  neglected  in  favor  of  terms  behaving  as  n  when  Eq.  (31)  was 

derived.  Equation  (31)  consistently  gives  slightly  too  high  a  growth  rate 

(i.e.  it  is  pessimistic),  a  feature  which  will  be  shown  to  persist  when 

variations  of  other  parameters  are  considered.  Notice  that  the  numerical 

solution  confirms  the  existence  of  a  peak  in  the  growth  rate  for  h  ■  1/yq>  as 

predicted  by  Eq.  (31),  again  independent  of  F,. 


The  significant  parameters  upon  which  the  small  signal  behavior  of  the 
layer  depends  include,  beside  the  externally  applied  bias  fields,  the  geometry 


factors  a,  b,  r^,  the  current  (v),  the  beam  energy  y  »  and  the  mode 

number  £.  We  proceed  to  consider  the  effect  of  each  of  these  separately  on 
the  growth  rate  of  the  longitudinal  mode,  w  •*  lwQ,  for  the  specific  case  of  an 
Astron-like  configuration,  h  ■  0. 

Unless  otherwise  stated  all  parameters  in  the  cases  considered  below  will 

take  the  following  nominal  "base  case"  values:  a  ■  0.5m,  b  -  2.2m,  r^  - 

0.99m,  r2  *  1.01m,  total  electrostatic  potential  difference  between  inner  and 

2  2 

outer  walls  ■  0.,  total  flux  between  inner  and  outer  walls  ■  ir(b  -  a  )  B 

oo 

-4 

where  BQO  -  48.2x10  T  is  the  field  required  to  hold  a  single  particle  at 

R  -  lm  with  y  -  3.  The  radial  density  profile  is  always  taken  to  be  parabolic 

and  symmetric  about  r«R»lm  with  specified  peak  value;  the  base  case  value  is 
7  -3 

5x10  /cc  which  gives  v/yq  "  3.94x10  .  The  base  case  azimuthal  mode  number 
l  is  1,  for  which  b+  +  b_  -  2.50  for  u  »  o)Q.  From  Eq.  (31)  the  normalized 
growth  rate  for  the  base  case  is  5.6%;  the  numerical  solution  of  the 
eigenvalue  problem  gives  5.35%. 

Figure  5  illustrates  the  comparison  of  growth  rates  as  calculated  by  Eq. 

(31)  and  by  a  numerical  solution  of  the  eigenvalue  problem  for  a  range  of 

1/2 

currents.  Over  the  wide  range  considered,  the  v  scaling  predicted  by  Eq. 
(31)  is  shown  to  hold  up  extremely  well  up  to  values  of  v/yQ  of  a  few 
percent.  Similar  excellent  agreement  is  generally  found  for  variations  in 
layer  thickness  (Fig.  6),  outer  wall  position  (Fig.  7),  inner  wall  position 
(Fig.  8),  particle  energy  (Fig.  9),  and  azimuthal  mode  number  (Fig.  10).  Some 
remarks  on  each  case  follow. 

In  the  case  of  varying  layer  thickness  with  the  maximum  density  fixed 

1/2 

(Fig.  6)  two  effects  are  competing;  these  are  the  basic  v  dependence  of  the 
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growth  rate  versus  the  stabilizing  effect  of  finite  thickness  (effectively, 
finite  frequency  spread).  The  finite  thickness  effect  is  clearly  second  order 
in  (t/R)  and  is  not  shown  explicitly  in  Eq.  (31)  but  it  _i£  included  in  all 
plotted  data  of  Figs.  5  -  10;  the  second  order  term  is  given  explicitly  in  the 
following  section.  In  the  examples  we  have  studied,  finite  thickness  effects 
have  been  small  and  have  not  been  effective  in  stabilizing  the  instability. 

In  Fig.  7  the  effect  of  outer  wall  position  is  illustrated.  As  the  wall 
is  moved  in  from  its  base  case  location  at  2.2m  the  growth  rate  is  observed  to 
increase  dramatically  for  a  while,  then  to  fall  off.  The  reason  for  the 
increase  is  the  approach  to  synchronism  of  particle  motion  with  a  cavity 
vacuum  mode,  that  is,  w  *  »  w  where  to  is  a  solution  to  b,  +  b  *  0; 

under  this  synchronous  condition  the  cavity  mode  fields  act  to  enhance  those 
established  by  the  dynamical  charge  bunching  due  to  the  negative  mass  effect. 
This  synchronous  case  has  been  given  its  own  name,  the  cyclotron  maser 
instability,  and  is  put  to  enormous  practical  use  in  the  gyrotron  family  of 
microwave  devices. ^  While  in  Figs.  5  -  10  we  have  consistently  evaluated 
b+  and  b_  at  to  ■  io)Q,  near  a  zero  of  b+  +  b_  this  is  clearly  inadequate  and 
Eq.  (31)  should  be  solved  as  a  cubic  polynomial.  Empirically  we  find  that 
evaluating  b+  +  b.  at  w  ■  is  adequate  when  b+  +  b_  >  l . 

As  the  outer  wall  is  moved  further  inward,  past  the  synchronous  point  the 
growth  rate  in  Fig.  7  drops  as  b+  +  b_  changes  sign.  This  drop  is 
attributable  to  a  "shorting  out"  of  the  azimuthal  field  E^  as  the  wall 
approaches  the  edge  of  the  layer.  An  identical  phenomenon  is  seen  as  the 
inner  wall  is  moved  outward,  Fig.  8.  Use  of  an  inductive  impedance  (b+  +  b_ 

<0)  to  stabilize  the  negative  mass  instability  has  been  proposed  by  Briggs  and 

32 

Neil.  We  remark  that  in  Fig.  8  no  synchronous  case  is  encountered  for  the 
parameters  considered,  as  the  inner  wall  is  moved. 


The  decrease  of  growth  rate  with  increasing  kinetic  energy  is  documented 
in  Fig.  9.  The  basic  reason  for  this  decrease  is  just  the  relativistic  mass 
increase:  the  azimuthal  response  of  a  particle  to  the  perturbed  field  E  is 
reduced  by  a  factor  yq^»  for  large  yQ.  No  synchronous  cases  are  encountered 
over  the  range  of  yQ  considered. 

In  Fig.  10  we  have  plotted  growth  rate  versus  azimuthal  mode  number  Z  . 
Agreement  is  good  between  the  dispersion  relation,  Eq.  (31),  and  the  numerical 
solution;  near  synchronous  cases  occur  for  Z  ■  5  and  9.  Though  predicted 
growth  rates  are  rather  large,  for  the  high  l  modes,  these  should  in  practice 
be  subject  to  stabilization  by  the  effects  of  finite  betatron  oscillations 
which  we  have  neglected  in  this  treatment.  In  any  event  we  expect  the 
dispersion  relation,  Eq.  (31),  to  begin  to  break  down  for  short  wavelengths, 
i.e.  £t/R  ~  0(1),  or  azimuthal  wavelengths  on  the  order  of  the  layer 
thickness.  A  WKB  treatment  of  Eq.  (14)  carried  out  for  this  case  yields  an 
eigenvalue  condition,  the  numerical  solution  to  which  would  appear  to  require 
more  effort  than  a  direct  numerical  solution  of  Eq .  (14)  itself. 

Finally,  we  report  that  in  the  initial  phase  of  our  investigation,  we 

attempted  to  derive  the  dispersion  relationship  by  constructing  a  quadratic 

form  for  the  differential  equation  (14),  under  the  assumption  that  <J>  ■  rE,  is 

approximately  continuous  across  the  E-layer.  The  latter  assumption  has  been 

9—1 1  17  18 

widely  used  and  has  been  considered  valid  ’  *  for  a  thin  relativistic 


electron  beam.  However,  this  line  of  investigation  led  to  an  incorrect 

2  2  > 

dispersion  relationship  if  Wp/oiQ  ~  0(1),  even  to  the  lowest  order  in  t/R.  In 
other  words,  to  account  correctly  for  the  DC  self  fields  in  the  present 


Eulerian  description,  the  tangential  AC  electric  field  should  not  be  assumed 


constant  across  the  E-layer,  regardless  of  the  thickness.  On  the  other  hand, 


our  dispersion  relationship  Eq.  (31)  correctly  accounts  for  the  self  fields, 


and  is  valid  for  arbitrary  beam  energy,  and  arbitrary  combination  of  EQ  and 
Bq,  as  already  stressed. 


IV.  Residual  Growth:  The  Dlocotron  Instability 

2 

Sufficiently  close  to  the  "zero  mass"  points  h  »  -  8Q/2,  h-*-  ±  », 

2 

Eq.  (31)  begins  to  be  dominated  by  the  neglected  terms  of  order  (t/R)  •  We 

2 

have  already  observed  this  phenomenon  near  the  point  h  =  -  B  /2  in  Fig  4.  In 

this  section  we  discuss  the  point  h  +  ±  »,  which  extreme  is  reached  in  the 

planar  limit:  r  -►  <*,  l  -*■  <*>,  l/r  =  k^  finite.  We  remark  that  the  vanishing  of 

the  negative  mass  growth  rate  in  the  planar  limit,  a  feature  of  Eq.  (31)  and 

clearly  expected  on  physical  grounds,  has  not  previously  been  demonstrated 

analytically,  to  the  authors'  knowledge. 

The  second  order  term  in  Eq.  (31)  is  given,  in  complete  generality, 

2  2 

by  (jJt/R)  mQh  where 


A’  ['2(b»  +  2h  +  (1+h)2) 
+2Cy?  (8?  +  2h)2  -  Y?(1+yV)(0?  +  2h) 2 


rj  (1  ±  h)  ,b+  ~  b~ , ) 2 
•■2  2, ,  2.  2.  '•b  +  b  ‘ 

V1  +  Y0h  )  + 


As  it  stands  A  includes  both  the  diocotron  instability  (£  terras)  and  the 
finite  thickness  stabilization  effect‘d  (E°  term)  referred  to  above.  The  last 
term  of  Eq.  (36)  contains  wall  boundary  effects;  specifically  we  may  recover 
the  stabilization  at  the  diocotron  mode  due  to  contact  of  the  layer  with 
perfectly  conducting  walls, ^  as  well  as  the  destabilizing  effects  of  finite 
wall  conductivity . 3A-37  we  aggurae  that  the  walls  are  many  wavelengths  away 
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1. 


a 

B 


W 


and  that  the  fields  are  electrostatic  then,  in  the  planar  limit,  = 

3  2 

Taking,  then,  the  planar  limit,  hioo-»-  -  eBQ/(mY0)  =  -  to  /y  ,  only  the  first 

2  L 

5  term  in  Eq.  (36)  survives  and  we  recover  the  classical  result4  for  the 

growth  rate  of  the  diocotron  mode 


Im(oj)  ■  |k  rv’/2| 
y  o 


(37) 


where  the  velocity  shear  vj  =»  .  Of  course  Eq.  (37)  may  be  obtained 

by  much  simpler  methods  than  those  employed  here;  our  point  is  only  that  it  is 
recoverable  from  the  present  formalism.  Note  that  the  dependence  on  the  line 


density,  v,  of  the  diocotron  growth  rate  is  just  v  ;  for  the  (non-synchronous) 

1/2 

negative  mass  instability  the  dependence  is  v  ;  for  the  synchronous  case  it 

.  1/3 

is  v 


The  relationship  between  the  diocotron  and  negative  mass  instabilities 

lO  1  O  OQ 

has  been  discussed  by  Neil  and  HeckrotteJO  and  by  Lau.  °  Mostrom  and  Jones 

have  recently  examined  the  electrostatic  case,  including  the  effects  of  shear 

29 

in  vz-  Davidson  and  Tsang  have  reported  analytical  and  numerical  results  in 
cylindrical  geometry. 


V.  The  Transverse  Mode 

An  electron  moving  in  a  field  satisfying  Eq.  (25)  where  Re(w)  is  a  vacuum 
guide  mode  may  be  shown  to  be  acted  upon  by  a  nearly  constant  electric  field, 
when  the  motion  is  averaged  over  its  gyro-orbit.  The  particle  therefore 
experiences  an  ?  x  $  drift  transversely,  toward  the  wall,  which  motion  brings 
the  particle  to  experience  yet  a  stronger  electric  field.  A  net  transfer  of 
energy  from  (to)  the  particles  may  be  shown  to  result  for  the  -  (+)  sign 
resonance  of  Eq.  (25).  This  mechanism  of  wave  growth  has  been  used  to  explain 
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the  operation  of  the  peniotron  oscillator. 


« 
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In  the  planar  limit  of  the  previous  section,  the  resonance  condition  for 
the  transverse  mode,  Eq.  (23),  becomes  y0(w-  kj,v0)  =  ±  ojc,  which  is  the  well- 
studied  mode  of  planar  magnetron  tube  theory. 1“7»33,34  xhe  factor  [Eq. 
(5)],  thus  again  appears,  as  it  did  in  Eq.  (31),  as  a  measure  of  the 
"planarity"  of  the  configuration. 

A  dispersion  relation  for  the  transverse  mode  may  also  be  obtained  from 
Eq.  (30).  Keeping  the  most  important  terms  we  have,  approximately 


|<b+  *  »->  -  5  TT  <lV  ‘»<K  +  ’>  ■  0  <38> 

o 


where,  near  the  zeroes  of  A 


-  (n2-  a2w2) 

a  ~  _ _ 

o  “  2  2 


1  n  2  UJ  .  2 

Z  (o  (1  -  8  - — ) 
p  ofcw 

o 


-(0 


q  = 


(1  +  h) 


Za  (1  -  02  -p-) 

o 


(39) 


(40) 


where  a  =  1  +  yq  h  .  Equation  (38),  using  Eqs.  (39)  and  (40),  agrees  with 
Eq.  (39)  of  Briggs  and  Neil^®  for  h  -  0,  to  leading  order  in  £ . 

Clearly  there  are  no  unstable  roots  of  Eq.  (38)  near  the  (simple)  zeroes 
of  Aq  unless  either  one  of  these  nearly  coincides  with  a  zero  of  b+  +  b_ 

(a  guide  mode  cutoff  frequency)  ££  b+  or  b_  contains  an  imaginary  part,  due  to 
finite  wall  resistivity  for  example  (see  Appendix  B).  Thus  the  transverse 
mode  (for  small  h)  depends  crucially  on  the  interaction  of  the  electrons  with 
their  external  surroundings,  unlike  the  negative  mass  and  diocotron 
instabilities,  the  mechanisms  of  which  operate  in  a  manner  that  Is  insensitive 
to  boundary  conditions  on  the  fields  at  distant  walls.  In  the  synchronous 


case,  with  perfectly  conducting  walls  the  predicted  scaling  of  the  growth  rate 

1/2 

of  the  transverse  mode  with  v  is  v  .  The  dependence  of  the  real  part  of  the 
frequency  of  the  mode  on  an  externally  applied  radial  DC  electric  field  may 
prove  useful  in  some  circumstances. 

The  peniotron  interaction  is  essentially  non-relativistic.  It  relies 
heavily  on  the  spatial  inhomogeneity  of  the  perturbed  wave  fields.  In 
contrast,  for  longitudinal  modes,  the  spatial  inhomogeneity  of  the  unperturbed 
motion  (i.e.  shear)  is  far  more  important.  Both  the  transverse  and 
longitudinal  modes  can  be  used  to  convert  the  rotational  energy  of  the 
electrons  to  rf  waves  efficiently,  however. 

In  the  planar  limit  Eq.  (38)  continues  to  predict  stability  in  the 
absence  of  a  resistive  wall  yet  it  is  well-known^--’  that  inclusion  of  a 
resonant  layer  satisfying  yQ(o)  -  k  vQ)  =  ±  in  the  beam  leads  to  wave 
growth.  The  resolution  of  this  contradiction  lies  in  the  failure  of  the 
Taylor  series  solution  to  Eq.  (14),  from  which  Eq.  (38)  was  obtained;  the 
resulting  growth  rate  is  non-analytic  in  %  and  one  mu3t  resort  to  numerical 
or  other  methods^  to  solve  the  eigenvalue  problem. 


VI.  Summary 

In  this  paper  we  have  attempted  a  general  treatment  of  the  linear 
stability  problem  for  laminar  electron  flow  In  cylindrical  geometry.  The 
basic  equilibrium  state  has  been  taken  to  be  maintained  by  radial  electric 
and/or  axial  magnetic  fields  (Gq.  (1))>  No  azimuthal  magnetic  field  nor  any 
axial  electron  motion  have  been  included  in  the  equilibrium  state.  The  linear 
stability  problem  for  azimuthal  and  radial  perturbations  has  been  formulated 
exactly,  fully  relativistically  and  fully  electromagnetically ,  including  all 
effects  of  self  fields.  The  stability  problem  reduces  to  an  eigenvalue 
problem  for  the  frequency  oi,  given  the  azimuthal  mode  number,  %  (Eq.  (14), 
with  associated  boundary  conditions). 

Our  efforts  have  been  focused  on  the  longitudinal  mode  (Eq.  (22)),  for 
thin  beams,  which  is  of  considerable  importance  in  accelerator  and  microwave 
device  research.  We  have  obtained,  and  favorably  compared  to  a  numerical 
solution  of  the  eigenvalue  problem,  a  dispersion  relation  in  the  thin  beam 
limit  (Eq.  (31))  which  applies  in  complete  generality  to  the  longitudinal  mode 
and  which  reproduces  all  classical  results  in  appropriate  limits.  Some 
interesting  differences  among  equilibria  regarding  the  negative  mass 
instability  have  been  pointed  out;  namely  we  have  found  a  simple  way  either  to 
stabilize  or  to  maximize  the  growth  of  this  mode.  This  finding  might  have 
practical  consequences  in  accelerator  or  microwave  tube  design. 

The  longitudinal  mode,  w  -  Iuq»Q,  encompasses  the  negative  mass,  electron 
cyclotron  maser,  and  diocotron  instabilities.  The  negative  mass  and  electron 
cyclotron  maser  effects  are  unique  to  cylindrical  geometry;  they  are 
fundamentally  relativistic  in  nature  when  the  motion  is  supported  solely  by  a 
magnetic  field.  They  are  even  more  pronounced,  especially  for  low  energy 
beams,  if  the  equilibrium  rotation  is  supported  solely  by  a  radial  electric 


field.  In  planar  geometry  both  of  these  instabilities  are  absent  and  only  the 

residual  diocotron  instability  remains  which  itself  may  be  stabilized  by 

placing  the  layer  in  contact  with  a  conducting  wall.^ 

The  transverse  mode,  td-£u)o  ■  ±  auQ  (a  is  defined  following  Eq.  (40))  has 

been  used  to  explain  the  operation  of  the  peniotron  device.  When  the  geometry 

is  cylindrical  (|y  h|«l)  this  mode  is  stable  unless  the  electron  motion  is 
o 

synchronous  with  a  cavity  mode  and/or  resistive  walls  are  present.  In  the 
planar  limit  aitiQ-*itsc/y  Q  and  we  identify  this  mode  as  the  Doppler  shifted 
cyclotron  resonance  considered  by  Buneman^*^  and  others^- »6>7»33  in  studies  of 
magnetron  operation.  This  mode  is  the  dominant  unstable  mode  for  planar,  high 
density  laminar  flow. 

Finally  we  remark  that  the  singularities  defined  by  Eq.  (24),  which  we 
have  not  examined  here,  may  be  worth  some  additional  study;  however  we  note 
that  in  both  the  vacuum  case,  5*0.  and  In  the  case  of  planar  Brillouin 
flow, 5, 33  ^2  m  J2'  the  singularity,  Eq.  (24),  is  removable. 
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Fig. 


Model  of  an  E-layer.  The  layer,  of  infinite  extent  in  z  (in  and 

out  of  the  page),  occupies  the  region  rj  <  r  <  r2  between  the  walls 

of  a  coaxial  guide  at  r-a  and  r"b.  The  electrons,  supported  by  an 

electric  field  E  (r)r  and  a  magnetic  field  B  (r)z,  move  in  concentric 
o  o 

circles  in  the  equilibrium  state,  either  clockwise  or 
counterclockwise  depending  on  the  equilibrium  type. 
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TWO  ROOTS 


corresponds  to  clockwise  or  counterclockwise  rotation  respectively  as 


A  ^ 

M 


Plot  of  the  dimensionless  "azimuthal  mass"  M  versus  the 
dimensionless  equilibrium  electric  field,  h.  The  actual  plot  shown 
is  for  the  case  y  *  1.5  but  all  axis  labels  are  expressed  generally 


in  terms  of  y  or  3 


numerical  solution  of  Eq.  (14);  the  dashed  line  Is  a  plot  of  Eq. 

_3 

(3t).  The  upper  pair  of  curves  is  for  v/yq  *  7.88  x  10  , 

_3 

C  *  1.42;  the  lower  pair  is  for  v/yq  *  1.57  x  10  ,  £  »  0.28. 
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ro/  (m)  III  I 


.14 


10 


Fig.  5.  Normalized  growth  rate  for  the  negative  mass  instability 

versus  v/yq  for  the  "base  case"  parameters:  a  *  0.5m,  b  -  2.2ra, 
r^  ■  0.99m,  ^  ■  1.01m,  yq  =  3,  l  -  1.  A  solid  curve  denotes  data 
obtained  from  a  numerical  solution  of  the  eigenvalue  problem;  a 
dashed  curve  denotes  data  from  the  dispersion  relation,  Eq.  (31). 
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Normalized  growth  rate  for  the  negative  mass  instability  versus 
layer  thickness.  All  parameters  take  their  base  case  values,  (see 
text)  except  that  the  peak  density  in  the  parabolic  profile 


ft 


4 


outer  wall  position.  All  parameters  take  their  base  case  values. 

(See  text.)  Near  synchronous  conditions,  Jtw  =oj  where  u  is  a 

o  v  v 

waveguide  mode  satisfying  b+  +  b_  =»  0,  the  negative  mass  instability 
is  strongly  enhanced  and  Eq.  (31)  should  be  solved  as  a  cubic 
polynomial.  The  synchronous  or  enhanced  negative  mass  instability  is 
often  called  the  cyclotron  maser  instability. 
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Normalized  growth  rate  for  the  negative  mass  instability  versus 
azimuthal  mode  number  £.  All  parameters  take  their  base  case 
values.  (See  text.)  A  +  denotes  a  result  from  Eq.  (31)  and 
an  x  denotes  a  result  from  a  numerical  solution  of  the  eigenvalue 


problem. 


Appendix  A 


Technique  for  Numerical  Solution  of 
Equilibrium  Problem 

In  this  Appendix  we  describe  a  simple  technique  based  on  Newton's  method 
which  we  have  used  to  solve  a  certain  eigenvalue  problem  associated  with  the 
calculation  of  laminar  E-layer  equilibria.  Given  the  electric  field  which  is 
trivially  solved  for,  having  specified  the  density  profile  and  potential  on 
one  wall,  we  must  find  the  momentum  profile,  u(r)  and  the  field  BQ(r)  subject 
to  the  constraint  of  specified  total  flux  Eq.  (7).  The  total  flux  might  be 
specified  in  an  experiment  in  which  a  beam  is  injected  and  contained  in  a 
chamber  for  less  than  a  magnetic  diffusion  time. 

Let  us  cast  the  problem  in  terms  of  the  dimensionless  fields,  and  a„ 

h,  15 

defined  in  the  text  following  Eq.  (6).  The  problem  then  is  to  find  u(r)  and 
Og(r)  subject  to 

b 

/  dr  Og  =  F  =  specified  constant  (A-l) 

a 


where  u  and  aD  satisfy 

15 


2  2 

l-£/y  -  y  h 

Ur  O  O  i 

u  "  rt — m — J 

2 

1  U) 

4  V'  u 

c 


(A-2) 


(A-3) 


2 

u 


,  2, ..1/2, 

OgCu  +1)  +  OgU 


(A-4) 


for  r^<r<r2*  Equation  (A-2)  guarantees  that  if  u(r^)  satisfies  Eq.  (A-4)  then 

u(r)  will  do  so  for  all  r.  The  algorithm  proceeds  as  follows:  An  initial 

guess  is  made  for  cig(r^)  using  the  value  of  the  externally  applied  B  field, 
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Appendix  B 

Evaluation  of  the  TE  Wave  Admittances 


In  this  Appendix  the  normalized  admittances  b^  referred  to  in  the  text  in 
Eqs.  (21a, b)  are  given  for  the  geometry  of  Fig.  1;  the  toroidal  and  planar 
cases  are  also  discussed. 


In  the  vacuum  regions  a<r<rj  and  r2<r<b  the  wave  equation  Eq.  (14)  is 


rdr*“TT - 2  if 5  + 


(  B— 1 ) 


<o  r 


-t 


the  general  solution  to  which  is 


♦  *  x  [ciJ*(x)  +  c2y^x>] 


(B-2) 


where  x  =■  wr/c,  and  are  the  standard  Bessel  functions,  and  C2  are 

constants,  and  in  this  Appendix  a  prime  will  denote  Note  that  the 

2  2 

"singularities"  at  x  =*  z  in  Eq.  (B-l)  are  only  apparent,  not  real,  i.e.  <)>  is 
analytic  at  these  points.  The  other  vacuum  fields  are 


f  -  iz 
lr  2  2  dr’ 

x  -Z 


Iz 


~r  E-\  • 

cZ  lr 


(B-3) 


(B-4) 


Using  the  definitions  in  the  text  Eqs.  (21a, b),  it  follows  that  b±  may  be 
generally  written 


Z  rCl+Jl(x2)  +  C2+Y£(x2) 
b  »  -  —  [  —  —  ] 

2  C1  JZ(X  2)  +  C2 vx2} 


(B-5a) 
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(B-5b) 


M 


b  -  i  [Cl  ^(Xl)  +  °2  ^(Xl>] 

X1  C1~J'(x1)  +  C2'Y'(Xl) 

where  2=  2^c* 

The  ratios  C^+/C2+  and  / C 2  are  determined  by  the  boundary  conditions 

at  r*b  and  r»a  respectively.  Some  special  cases  of  interest  are: 

1.  Perfectly  conducting  wall  at  r»rw:  (rw=*  a  or  b) 

ci/c2  ■  -T;  <v  <B-6> 


2.  Wall  with  (complex)  dielectric  e(a»): 


Cl/C2 .  -  ,yv  -  (W) 

J£<V  "  (?/co)Ji(xw)d7lnZ(y) 

1/2 

where  5  is  the  surface  Impedance,  (y/e)  »  CQ  “  376.70, 

1/2 

y”  <r  ,  (ey)  ,  Z(y)  -  J  (y)  for  inner  wall, 

w  JC 

H^^(y)  for  outer  wall  where  is  the  Hankel 

function,  e  and  y  refer  to  the  wall  material. 


3.  No  inner  wall: 


C2  “° 


A.  No  outer  wall: 


(B-8) 


Cl+/C2+  "  -  1 


(B-9) 
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5.  Electrostatic  limit;  perfectly  conducting  walls  at  r-a,b: 


B  +  1  B  +  1 

b+  -  3-^  ;  b.  -  3^  (B-10) 

where  B+  =  (b/r2)2^  and  B_  =  (r^/a)2*.  Note  that  b^.  >  0. 

6.  Planar  limit  (  r-n»,  £-m»;  j,/r=ky,  b-a,  b-r^,  b-^,  r2"ri»  al* 
remain  finite);  perfectly  conducting  walls: 

b+  -  -  (k  la)  cot  [a(r1~a) ]  (B-ll) 

b_  -  -  (ky/a)  cot  [a(b-r2)]  (B-12) 

where  a  =  (w2/c2-  k2)1/2.  The  planar  limit  of  (B-10)  is 
just  the  electrostatic  limit  of  (B-ll)  and  (B-12),  as 
expected,  in  which  case 

b+  -  coth  [ky(b-r2)J  (B-13) 


b_  *  coth  [ky(r^-a)J.  (B-14) 

When  a  wall  is  resistive  the  resulting  dissipation  is  represented  by  an 
imaginary  part  in  the  corresponding  admittance  b±.  (In  Eq.  (B-7)  e(w)">i o/w 
for  a  good  conductor  of  conductivity  a.)  That  such  dissipation  can  lead  to 
disruptive  beam  instabilities,  even  for  a  "positive  mass”  beam,  has  been  known 

1  e 

since  the  pioneering  work  of  Neil  and  Sessler  and  Laslett,  Neil,  and 
Sessler.^  It  is  in  fact  these  resistive  modes,  rather  than  the 
(fundamentally  dynamical)  negative  mass  instability  which  are  thought 
ultimately  to  limit  the  beam  current  in  cyclic  electron  accelerators. 
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If  a  wall  Is  not  smooth  but  contains  some  structures  (cavities,  fins, 

etc.)  a  common  practice  is  to  calculate  an  approximate  value  for  the 

admittance  and  to  use  it  as  a  boundary  condition  for  some  approximation  to  Eq. 

(14)  in  the  text.  In  this  case  however  the  problem  is  not  being -treated  fully 

self  consistently  since  the  equilibrium  of  Section  II  would  not  be  strictly 

correct,  i.e.,  the  correct  equilibrium  would  no  longer  have  azimuthal  symmetry 

(and  would  be  much  harder  to  calculate). 

For  the  case  of  a  perfectly  conducting  wall  we  have  plotted  in  Figs. 

B-1,2,3,4  the  quantity  b+  +  b_  evaluated  for  x^  ■  X2  ■  for  various  values 
2 

of  E^  =  me  (yo~1),  assuming  a  thin  layer.  It  is  this  quantity,  b++  b_,  which 
enters  the  dispersion  relation  for  the  longitudinal  mode  discussed  in  the 
text. 

For  a  toroidal  (accelerator)  geometry,  the  dispersion  relation  Eq.  (31) 
is  expected  to  be  replaced  by 

nj  -  -  eog*2  (2v/Yo)(c2/R2)(mc2/a>o)Omo/3E)ext  (B-15) 

where  (3u>0/3e)ext  denotes  the  derivative  of  the  circulation  frequency  with 
respect  to  total  particle  energy  evaluated  as  if  the  particle  were  acted  upon 
only  by  the  external  electric  and  magnetic  fields.  In  Eq.  (B-15),  we  have 
identified  the  geometric  factor  g  of  refs.  [8,11]  with  l/(fceo(b+  +  b_)) ,  as 
pointed  out  in  ref.  10.  Note  that  the  dispersion  relationship  (B-15)  includes 
self  field  effects  and  that  g  is  always  positive  for  toroidal  geometry  (with 
smooth,  perfectly  conducting  walls). 


» 


1 

I 


12  3  4 


Fig.  B-l.  Normalized  admittances  vs.  outer  wall  position  for  various  inner 

2 

wall  positions  for  &  -  1  and  me  Cto~1 )  ■  1  keV. 
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b+— >  ±  oo 


Fig.  B-3.  As  in  Fig.  B-l  but  for  me  (y  -1) 

o 


40  keV. 


IT -0.2  (b_*  1,32) 

R  ~ 

2 

As  in  Fig.  B-l  but  for  me  (,y  -\)  »  300  keV.  The  zeroes  of  b+  +  b. 
correspond  to  waveguide  cutoff  frequencies. 
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